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Abstract 

The NLO term in the OPE of the quark propagator vector part and the 
vertex function g\ of the vector current in the Landau gauge should be dominated 
by the same (A 2 ) condensate as in the gluon propagator. On the other hand, the 
perturbative part has been calculated to a very high precision thanks to Chetyrkin 
and collaborators. We test this on the lattice, with both Wilson-clover and GW 
(Ginsparg- Wilson) overlap fermion actions at (3 = 6.0, 6.4, 6.6, 6.8. Elucidation of 
discretisation artefacts appears to be absolutely crucial. First hypercubic artefacts 
are eliminated by a powerful method. Other, very large, non perturbative, 0(4) 
symmetric artefacts, impede in general the analysis. However, in two special cases 
with overlap action - 1) for Z$ ; 2) for g±, but only at large p 2 - we are able to identify 
the (A 2 ) condensate; it agrees with the one resulting from gluonic Green functions. 
We conclude that the OPE analysis of quark and gluon Green function has reached a 
quite consistent status, and that the power corrections have been correctly identified. 

1 Laboratoire associc au CNRS, UMR 8627 



A practical consequence of the whole analysis is that the renormalisation constant Z^ 
(= Z^ 1 of the MOM scheme) may differ sizeably from the one given by democratic 
selection methods. More generally, the values of the renormalisation constants may be 
seriously affected by the differences in the treatment of the various types of artefacts, 
and by the subtraction of power corrections. 



PACS: 12.38.Gc (Lattice QCD calculations) 
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1 Introduction 



The study of the quark propagator and vertex functions in momentum space has been extensively 
pursued in the literature starting in the 70's with analytical considerations [1-3], later completed by 
the discovery of the presence of a contribution of the A 2 operator, due to gauge fixing ( [4]). Numerical 
lattice QCD has more recently treated this issue [5-12]. The scalar part of the quark propagator is 
related via axial Ward identities to the pseudoscalar vertex function. The role of the Goldstone boson 
pole in the latter has been thoroughly discussed [13-16]. 

In this paper we will mainly concentrate on the vector part of the inverse quark propagator, the 
one which is proportional to Z^, and on the vector vertex function gi related to it by the Ward 

identity, or, equivalently, on the ratio Zy OM (p) = ^jpj ■ One of the reason is to check for the effect 

of the (A 2 ) condensate which has been discovered via power corrections to the gluon propagator 
and three point Green functions at large momenta [17-20] (let us recall that, in order to perform 
trustable perturbative calculations, we take always p > 2.6 GeV). OPE shows that the effect of the 
(A 2 ) condensate should have almost the same magnitude in the Z^(p 2 ) of the quark as in the gluon 
propagator. Moreover, the perturbative-QCD corrections are known to be varying especially slowly 
(the anomalous dimension being zero at one loop in the Landau gauge), as we shall recall later. This 
should give a favorable situation to display the power corrections (recall that in the gluon case, it was 
difficult to disentangle the power and the logarithmic corrections, which were moreover very sensitive to 
the value of Aqcd)- Now, for the Wilson-SW action, the crude values plotted in litterature (with only 
a selection of democratic points) for Z^{p 2 ) are extremely flat above 2 GeV [24] This was formerly 
considered as being natural, as a consequence of the vanishing of one-loop anomalous dimension. 
However, thinking more about it, it must be considered on the contrary as very worrying since it 
means the decrease predicted both by the perturbative QCD corrections and the (A 2 ) condensate 
would not be seen. This lead us to start a very systematic study of the problem, with the following 
series of improvements on earlier works: 

• We reach an energy of 10 GeV by matching several lattice spacings so we are in a better 
position to eliminate lattice artefacts and also to identify power corrections, which requires a 
large momentum range. 

• We make use of a very efficient way of eliminating hypercubic artefacts which we have first 
elaborated while studying gluon propagators [21], [22]. In a recent paper [23], this method 
has still been improved for the specific case of the quark Green functions, where such artefacts 
are huge, especially in Z$ and for the overlap action. Hypercubic artefacts have often been 
cured by the "democratic" method which considers only momenta with equilibrated values of 
the components. This method points in the right direction but, as we shall discuss in more 
details, is by far insufficient when the hypercubic artefacts reach such a level as illustrated in 
this article. 

• We make a systematic use of Ward-Takahashi identities relating the quark propagator and the 
vertex function. According to these Ward identities, Zy OM should be independent of p 2 up to 
artefacts. Then, Zy OM is a very sensitive test of the presence of artefacts. The fact that we 
observe a strong dependence of the lattice Zy OM on p 2 , shows unambiguously the existence 
of large remaining discretisation artefacts, this time respecting 0(4) symmetry, and decreasing 
with negative power of momentum. We can also check the consistency of Zy OM with other 
determinations of Zy . 
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• The above mentioned 0(4) symmetric artefacts could constitute a very serious problem, since 
we do not know how to eliminate them, and they would hide the OPE power corrections. For- 
tunately, in the case of overlap fermions (with s = 0), we find that for Zj,, there remains only 
very small 0(4) artefacts. Then, we can obtain a satisfactory estimate of the 1/p 2 correc- 
tion due to the (A 2 ) condensate. For this reason, we consider mainly overlap fermions (with 
s = 0), although they have some specific unconveniencies (for example very large perturbative 
corrections). 

This work, as well as the preceding one, clearly shows that lattice artefacts are overwhelming at the 
start, both hypercubic and 0(4) ones. Hypercubic ones have been shown to be cleanly eliminated by 
our method. Remaining 0(4) ones, on the other hand, cannot, and we do not foresee the possibility of 
a similarly efficient method, wherefrom we have to rely on situations where they are small for reasons 
which are not known, so that their smallness appears accidental. In that respect, it may seem worrying 
; but, on the other hand, we are very happy to have found that the OPE can be checked to a good 
accuracy. 

Another obvious interest of the study is then to improve the determination of the MOM renor- 
malisation constants, by taking into account both continuum power corrections and artefacts ; and, 
indeed it indicates that much care must be exerted in using MOM renormalisation approach when high 
precision will be required, a point which has already been illustrated by the Goldstone contribution to 

gMOM 

In section El we will recall some theoretical premises, in section El we will indicate the lattice 
conventions and the simulations which we have performed, in section [4] we will recall briefly the 
method to eliminate hypercubic lattice artefacts, in sectionsOOEl we discuss other artefacts (Lorentz 
scalar artefacts, volume effects), in sections 171 151 191101 we will give the results and, in section 1TT1 we 
will give our conclusions and further discussions. 

2 Theoretical premises 

We work in the Landau gauge. Let us first fix the notations that we will use. We will use all along the 
Euclidean metric. The continuum quark propagator is a 12 x 12 matrix S(pn) for 3-color and 4-spinor 
indices. The inverse propagator is expanded according to : 

S- 1 (p) = 5 a ^(p 2 ){ii> + m(p 2 )) (1) 

where a, b are the color indices. Z^{p 2 ) being a standard lattice notation (for the precise lattice 
definition, see below, section |3J). Obviously, one has in the continuum, with trace on spin and color: 

Z^p 2 ) = 1/12 Tr(S(ri lflPf ,)/p 2 (2) 

Sometimes, one uses the alternative quantity: 

b(p 2 ) =Z lp (p 2 )m(p 2 ) (3) 

to describe the scalar part of the propagator. 

Let us consider a colorless vector current qj^q. The three point Green function is defined by 

C>, q) = J d 4 xd 4 y e^+^' < q(y)q(x) lfl q(x)q(0) > (4) 
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The vertex function is then defined by 

r M (p, q) = S-\p) G^{p, q) S^ip + q) (5) 

In the whole paper, we will restrict ourselves to the case where the vector current carries a vanishing 
momentum transfer q^. In the following we will omit to write q^ = and we will understand 
as the bare vertex function computed on the lattice. 
From Lorentz covariance and discrete symmetries 

r^(p) = s a ,b [gi{p 2 hn + m{p 2 )p^ + gs(p 2 )p^ + m(p 2 )h»,j>)\ (6) 

which should be obeyed approximately on the lattice, as we checked. 



2.1 Renormalisation and Ward-Takahashi Identities 

The renormalised vertex function is then Z^ZyY '^(p). Here, we must say something of conventions for 
renormalisation constants. The standard definition of renormalisation constants has been to dividethe 
bare quantity by the renormalisation constant to obtain the renormalised quantity (except for photon 
or gluon vertex renormalisations Z\ which we do not use). Z 2 is the standard renormalisation of 
fermions ipbare = V~Z~2ipR, S(p) = Z2([J 2 )Sn(p). In principle renormalisation of composite operators, 
for instance Zy, should be defined similarly. We have followed this convention in our works on gluon 
fields. But, in the case of quark operators, an opposite convention has become standard in lattice 
calculations : (^pOip)b are = Zq (^}0^)r ; we feel compelled to maintain this convention for the 
sake of comparison with parallel works on the lattice. This explains our writing of the renormalised 
vertex function. In the continuum Zy = 1 (conserved current). We keep Zy since the local vector 
current on the lattice is not conserved, and the discrepancy, which is of course an artefact, generates 
however finite effects in graphs due to additional divergencies multiplying the a terms (which have 
higher dimension). 

The Ward identity in the renormalised form tells us that at infinite cutoff : 



1 ■ -i^K * (') 



After multiplying both sides by Z 2 1 to return to bare quantities 



Zvr»(p) = -i-^s- 1 (p), (8) 



which from Q-© implies 



Z vgi (p 2 ) = Z^p 2 ), Zyg 3 (p 2 )=2^Lz^p 2 ), -Zyg 2 (p 2 )=2^Lb(p 2 ), g 4 (p 2 ) = (9) 

We note that the first equation El implies that Zy is independent of the renormalisation scheme up 
to artefacts(\t is a ratio of bare quantities). Of course, this will hold up to terms vanishing as inverse 
powers of the cutoff at infinite cutoff, which are called artefacts in the lattice language. It must be 
recalled that, on the lattice, the Ward identity is not exact, but holds only up to artefacts, because 
we work at finite cutoff, and the deviation will be found very large in some cases. A very important 
consequence of the Ward identity for our study is that the ratio gi(p 2 )/Z^p(p 2 ) is constant up to 
artefacts, or that deviations of this ratio from a constant are pure artefacts. 
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Defining analogously the vertex function of the pseudoscalar {qj^q) density, 



r 5 (?v) =g 5 (p 2 h 5 (io) 

the axial Ward identity implies 



^rn q g 5 (p 2 ) = b(p 2 ) (11) 



where m q is 1/(2 a) (1/k — 1/k c ). 



2.2 MOM renormalisation ; radiative corrections 

To perform renormalisation on the lattice, we appeal as usual to the convenient MOM schemes, which 
does not refer to a specific regularisation. To speak technically, the precise renormalisation scheme 
that we use is the one called Rl' by Chetyrkin, eq (26) in ref. [25]. This is in fact the most standard 
MOM scheme in the continuum, developed a long time ago by Georgi, Politzer, and Weinberg. It 
consists in setting the renormalised Green functions to their tree approximation at the renormalisation 
point p 2 = /j, 2 , in the chiral limit. The inverse bare propagator S^ 1 (p) is normalised through : 

S^(p)\ p 2 = ^ = 5 a ^ b {i^ + m(p 2 )) (12) 

Making p 2 = /j, 2 shows that : 

Z 2 MOM (p 2 ) = Z^p 2 )^ (13) 

up to artefacts. To renormalise the bare vertex function gi(p 2 ), we multiply it by the factor 
Z^ OM \n 2 )Zy OM (/j, 2 ). The MOM renormalisation for g 1 must be chosen so that the renormalised W- 
T identity (eq. U$ gf (p 2 ) = Z^(p 2 ) holds, therefore gf (p 2 = pi 2 ) = 1 ; we deduce that Z^ OM {pL 2 ) is 
the ratio Z^(/j, 2 ) / g\(n 2 ) up to artefacts, but as remarked above, this must be nothing but the scheme 
independent Zy : therefore this ratio Zy OM (/j, 2 ) is independent of fi 2 , up to artefacts. From now 
on, we define Zy as the ratio Z t p(p 2 )/gi(p 2 ) (measured in fact on the lattice); we write : 

Z^ OM {a~\p 2 ) = Z^p 2 )/ gi (p 2 ) (14) 

which recalls that Zy , which should be in fact independent of p 2 in the limit where the cutoff 
a -1 is infinite, is not so at finite a -i.e. there are artefacts. 

The Ward identity implies that Z^(p) and g\(p) have in particular the same perturbative p scale 
dependence. From the calculations of Chetyrkin et al. [25] [26], we may express, for example, the 
perturbative running of Z^ at large p as a function of the running omom(p)- This is our main choice 
throughout this paper, although we discuss the effect of substituting an expansion in asjjg(p). More 
precisely, we will always choose to use the definition of «mom(p) by the symmetric three-gluon vertex. 
The advantage of quarks is that we can reach an accuracy of four loops in the RG expansion, because 
we do not need j3^ IOM (for symmetric «mom(p)), since the dimension of the fermion is at lowest 
order in the Landau gauge. Since the expression is lengthy and not necessary for present understanding, 
we refer the reader to the appendix|Al Of course, even with such accuracy, such an expression cannot 
be expected to hold for too small p : we esteem the lower bound to be p m in = 2.6 GeV, from our 
experience in the case of gluons ; indeed, we must avoid to go down too much close to the "bumps" 
which manifest clearly that the gluon Green functions become non perturbative. The perturbative 
calculation requires a value of Aqcd i one advantage of Z^(p) 2 is that it is not so much sensitive to 
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this value as the gluon quantities, because of the vanishing of the LO fermion anomalous dimension 
in the Landau gauge ; we choose A^r = 0.237 GeV from a previous analysis ( [17] ), not far from 
the ALPHA estimate [28], = 0.238(19) GeV ; we will discuss in the end the sensitivity of our 
results to this choice. 

Z^(p 2 ) and gi(p 2 ) should have also the same non perturbative power corrections, up to a constant. 
We consider them now. 



2.3 Power correction from the < A 2 > condensate 

An OPE analysis as those performed in refs. [17-20] leads to consider a (A 2 ) condensate coupled to 
the quark propagator and vertex in Landau gauge. Let us recall that such a condensate could not 
contribute to gauge invariant Green functions, and is present only in (gauge fixed) gauge non invariant 
Green functions. The meaning and magnitude of such a condensate has been extensively discussed in 
the recent literature. Our aim here is to detect its effect on the Green functions through OPE, which 
provides a way of testing theoretical ideas on its existence and magnitude. 
For the propagator, we can write: 



S-HP) = S£t(p) + ^i> pt^l ^ + ■ ■ ( 15 ) 

where we only keep the leading term in The calculation of the coefficients of the OPE has been 
performed in the chiral limit, and therefore one has as far as possible to stay near this limit. 

In the renormalisation prescription denoted by " Rl' " in Chetyrkin papers (which amounts to 
the standard MOM of Georgi and Politzer in the chiral limit), and expanding everywhere in terms of 
af OM , we obtain from Eq. flEJ) (see appendix!!]): 



(0) 

Z^(P 2 ) gTV) 32^ /oM\-V < (A 2 )M > 1 
Z^V) S^V) 3 WaOJ 4(JVg-l) p 2 [ ' 

where S^ 1 = ip* 5 a ^ Z^ and S~,} Tt = itf 5 a ^ Z^ CTt up to 0(m q ^>/p 2 )-terms. The condensate < 

(A 2 )r([i) > is renormalised at the scale \i. Z^ 7t is given in eq. (J21) and the coefficients (70 being the 
fermion anomalous dimension to lowest order) are 

a 11 n (°) 35Nc 35 f~i t\ 

/?o = 11, 70 = Y A 2 = —^j- = -j (17) 

As we have noted, Zy OM (p 2 ) should be constant in p 2 from the Ward identity ® , up to artefacts 
; then it cannot receive any power correction from A 2 , and therefore, g\ receives exactly the same 
contribution from the condensate as Z^. We will use this as a very useful test. 

The essential step is then to fit this formula on the lattice data to extract (A 2 ). The renormalisation 
constants at each j3, Z^{^l 2 ,I3) will enter in the fit as free parameters to be determined, although 
they would be expected a priori to be close to lattice perturbation theory predictions. Of course, in 
general, we have to add lattice artefacts to ea.([16|), and one of the main problems we will discuss is 
how to determine them accurately. 

An important warning must be made here, concerning the low accuracy in the perturbative calcu- 
lation of the Wilson coefficient of A 2 written above : namely, it is only tree order with renormalisation 
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group improvement. Expanding in terms of af ° , although it may seem natural, is completely 
arbitrary, and one would wish the results to be the same with a^ s . While this is the case to a good 
precision for Z^ CTt , this is obviously not the case here, due to the low order of the expansion : a s (p) is 
quite different in the two schemes : at p = 2.6 GeV, the ratio MOM /MS is around 2 and decreases 
slowly down to 1.4 at 10 GeV ; taking into account the anomalous dimension amounts roughly to 
replace a s (p) by a s (/j,) ; then the ratio of coefficients in terms of the two coupling constants is only 
slightly closer to 1 : it is 1.5 in average over the whole range . This means that the coefficient is 
reduced by 50% when using a^f s . This is due to the fact that ratio of coupling constants decreases 
only very slowly up to the largest available momenta. As a consequence, the determination of (A 2 ) 
obtained by fitting the lattice data will be automatically affected by the same amount. We will give 
the results with the convention of using everywhere a^ 10M , as we have done for gluons. We shall first 
show that the power correction is indeed present and well determined, and then express it in terms 
of the condensate value, which suffers from the above uncertainty. We also note that the ratio of 
condensates fitted from gluons and quark Green functions, which should be 1 ideally, is not affected 
by this uncertainty, since the Wilson coefficients relative to the various Green functions differ mainly 
by purely algebraic numbers (the anomalous dimensions differ only slightly) 1 . 



3 Lattice calculations 

We have first used SW-improved Wilson quarks (often called clover) with the csw coefficients com- 
puted in [32]. 100 quenched gauge configurations have been computed at (3 = 6.0,6.4,6.6,6.8 with 
volumes 24 4 , 16 4 and 8 4 . We have performed the calculation for five quark masses but in practice, 
for what is our concern in this paper, the quark mass dependence has not surprisingly proven to be 
negligible ; anyway, since the theoretical calculations are performed in the chiral limit, we have to 
work as close as possible to the chiral limit ; then, we present only for simplicity the results for the 
lightest quark mass, about 30 MeV, i.e. 

k = 0.1346, 0.13538, 0.13515, 0.13489 for = 6.0, 6.4, 6.6, 6.8 (18) 

It should also be mentioned that all the results presented for clover action refer to the 24 4 lattices 
unless stated otherwise. 

In addition to improved Wilson fermions, the use of overlap fermions [33] has revealed necessary, 
and even crucial to obtain a good determination of the power correction. We have used approximately 
the same physical masses i.e. as in the improved Wilson case 

am = 0.03, 0.01667, 0.01. for = 6.0, 6.4, 6.6, 6.8 (19) 

with s = and volumes of 16 4 . The bare mass itiq and s are defined from 

Awcr = (1 + s + am /2) + (l + s- am /2) ^_^ w( ~ (1 + s)) == === (20) 

v /,D w (-(l + s))tZ) w (-(l + s )) 

where D w (— (1 + s)) is the Wilson-Dirac operator with a (negative) mass term — 1 — s 

Av(-1 -s) = ^(V M + V*) - ^aV*V„ - 1 - s , (21) 



1 While finishing the article, we have become aware, thanks to D. Becirevic, of the calculation of the 
two-loop anomalous dimension of A 2 by [27], in the MS scheme 
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We use s = 0. s = 0.4 is considered preferable from locality requirements [34], however the 
difference is slight as soon as j3 is larger than 6.0. The reason for using only the small lattice 16 
is well known ; it is due to limitation in the special treatment needed for small eigenvalues of the 
Neuberger operator. In practice, as for clover action, we discuss only the lighest quark mass, roughly 
corresponding to the same 30 MeV. 

The propagators S(x, 0) from the origin to point x have been computed and their Fourier transform 

S{p)=XV *•*£(*,()) (22) 

X 

have been averaged among all configurations and all momenta p^ within one orbit of the hypercubic 
symmetry group of the lattice, exactly as for gluon Green functions in [21]. 

In the case of overlap quarks the propagator and other Green functions are improved according 
to a standard and exact procedure [35] which should eliminate 0{a) discretization errors in Green 
functions, at large p, in the perturbative regime 2 : 

S.(p) = f iP) - % (23) 
1 — amo/2 

From now on, the notation Sip) will represent the improved quark propagator in the case of overlap 
quarks and the standard one in the case of clover quarks. 
In both cases we fit the inverse quark propagator by 

S-\p) = 5 a , b Z^ (p) {i~j> + m(p 2 )) (24) 

according to eq. (JJ) and where p^ is defined in eq. (|28|). We write Z^(p) because of the loss of the 
Lorentz invariance. 

Z^(p) can then be written as : 

Z 1 p(p) = ^Tv[ lll p,S- l (p)}/(p) 2 (25) 

The three point Green functions with vanishing momentum transfer are computed by averaging anal- 
ogously over the thermalised configurations and the points in each orbit 

G>(P, Q = 0) =< 75 S(p)^ 75 7^ S(p) > . (26) 

where the identity S(0,x) = 7s5^(x,0)75 has been used. The vertex function is then computed 
according to eq. © and we choose for the lattice form factor gi : 



91(P 2 ) = 



F^(p,q = 0) ( 7 M 



(27) 



where the trace is understood over both color and Dirac indices. 

Finally, according to the Ward identity we compute Zy OM simply from eq. (|14|) where any 
effective p 2 dependence of Zy should come only from lattice artefacts. 

Throughout this paper we will use the values in the following table[Jfor the lattice spacings, which 
follow the j3 dependence found in ref. [36], appendix C, formule C.l, 

2 Indeed, in our opinion, the argument on the vanishing of artefacts uses chiral symmetry of vacuum 
matrix elements, which holds only when spontaneous symmetry breaking can be neglected. 
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p 


6.0 


6.4 


6.6 


6.8 


a- 1 (GeV) 


1.966 


3.66 


4.744 


6.1 


a (fm) 


0.101 


0.055 


0.042 


0.033 



Table 1: Lattices spacings 

4 Elimination of hypercubic lattice artefacts 

4.1 Classification of artefacts 

The question of eliminating lattice artefacts has been perhaps our main difficulty in this work. Since 
we will have a detailed discussion, it is useful first to remind the main species of artefacts which are 
expected. First, we have discretisation artefacts 3 , which themselves split into two : i) hypercubic 
artefacts, which are the most visible because they break the elementary 50(4) symmetry ; they are 
seen, as we plot invariants of Green functions as function of p 2 , as a large discrepancy between the 
value for different orbits at the same p 2 . With some simple treatment, it is easy to get one relatively 
regular function of p 2 . Nevertheless, in general, there remain non analytic oscillations, and to eliminate 
them is both important to obtain the final physical result, and demanding sophisticated methods. 

ii) £0(4) invariant discretisation artefacts, which remain after elimination of the cubic ones, and 
which will be discussed later. Let us say that this is the weakest point, because we do not have 
theoretical principles to determine their form, neither there is a systematical empirical method to 
determine them. 

Still, there may be finite volume artefacts which will be discussed also later, very shortly since they 
do not seem sizable. 

4.2 Hypercubic artefacts 
4.2.1 Generalities 

In successive papers [37,38], we have elaborated a very powerful method to deal with hypercubic 
artefacts i.e. with those discretisation artefacts which come from the difference between the hypercubic 
geometry of the lattice and the fully hyper-spherically symmetric one of the continuum Euclidean space. 
The principle of this method 4 is based on identifying the artefacts which are invariant for the H± 
symmetry of the hypercube, but not for the SO (4) symmetry of the continuum. 

Let us set the problem more precisely. Since we use hypercubic lattices our results are invariant 
under a discrete symmetry group, H4, a subgroup of the continuum Euclidean SO (4), but not under 
SO(4) itself. This implies that lattice data for momenta which are not related by an H4 transformation 
but are by a 50(4) rotation will in principle differ. Of course this difference must vanish when a —> 
but it must be considered among the discretisation effects, i.e. ultraviolet artefacts. For example, in 
perturbative lattice calculations one encounters the expressions 

= - sin {~y) , Pv = \ sin • ( 28 ) 

3 We will use the term "discretisation artefact" preferably to the other common one, "ultraviolet arte- 
fact", because, as we shall find, these artefacts may show up at small p as well, due to non perturbative 
effects. 

4 The initial idea is due to Claude Roicsncl 



9 



Both are equal to p^ up to lattice artefacts: 

p 2 ^T,pl=p 2 - T2 a2p[A] + • • • ' p2=p2 ~ r 2pW + • • • ' where ^ [2n] - E pl n (») 

p 2 , p 2 , p 2 , a 2 p^ are invariant under H4 but only p 2 is under SO (4). For example, the momenta 
2x/L(l, 1, 1, 1) and 2tt/L(2, 0, 0, 0) have the same p 2 but different pM ,p 2 and p 2 . In other words, 
if we call an orbit the set of momenta related by H4 transformations, different orbits correspond to 
the same p 2 . In general different orbits have different pW. The hypercubic artefacts can be detected, 
considering a given quantity at a given p 2 , by looking carefully how it depends on the orbit. 

These hypercubic artefacts, sometimes called "anisotropy artefacts", have been a long standing 
problem in lattice calculations, and of course, methods have been devised since a long time to handle 
them. The general idea has been the so-called "democratic" one : the hypercubic effects are minimal 
when the four components of the momentum for a given p 2 do not differ too much (this is democracy 
between components ; ideal democracy is for diagonal p oc (1,1,1,1)). Then the question is how to 
make this criterion quantitative, the rationale being to find a compromise between two contradictory 
requirements : 1) to be as much democratic as possible, which tends to reduce the number of points 
2) to retain enough points to have a real curve. 

The precise criterion is often something of a secret recipe, not communicated in papers. On the 
other hand , studying the gluon propagator, the authors of ref. [39] have made explicit a selection 
method, keeping only the orbits having a point within a cylinder around the diagonal. Several other 
similar criteria have been written. 



4.2.2 Our method : the extrapolation method 

The alternative idea which we proposed, on the contrary, relies on the use of all the orbits, and a 
method to extract the physical point from an extrapolation of the different orbits. A first successful 
application of this idea was for the gluon propagator [37]. We recall here the final refined form of the 
method, the so-called pl 2n ^ extrapolation method, presented in [38], and which has been shown to be 
necessary to obtain satisfactory results for Z^. 

In order to perform a global fit we start from the remark that in this paper we are dealing with 
dimensionless quantities, g\ and Z^. It is thus natural to expect that hypercubic artefacts contribute 
via dimensionless quantities times a constant 5 . Next we assume that there is a regular continuum 
limit. We denote the generic Green function as Q, which depends a priori on p 2 , but also a 2 p^\ 
a A p^ : Q(p 2 , a 2 p^\ a 4 p^\ a 4 p' 4 l, • • • ), and we Taylor expand it around the 0(4) symmetric limit 
Q(p 2 , 0, 0, 0, • • • ). Of course we must truncate this Taylor expansion of Q in a and we choose to 
expand it up to a 4 . Note that at this stage, the function Q(p 2 , 0, 0, 0, • • • ) may still depend on a 
through terms of the form a 2 p 2 , etc., but we do not consider presently this further dependence. 

The lattice results are H4 invariant and thus typically functions of p 2 , p 2 and p-p 6 . Then, let 
us consider a typical H4 invariant and dimensionless term and expand it using eq. J2U): 



2k (p 2 ) l (p 2 ) m (p-pr 
(p 2 y (p 2 ) m ' (p ■ p) n> 



a 



(«v) fc [i+«i/ 7 + ... n « 2niP —^+ -i ( 3 °) 



5 We neglect a possible logarithmic dependence on p 2 . 
6 In all this discussion we consider the mass as negligible. 
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where I + m + n — (I' + m' + n') = k. In order to have a continuum limit k > 0. As we expand 
the serie up to a 4 we have k + n j ^ 2. So for k = we obtain the following artefacts: a 2 ^- , 

a 4 (^2-) 2 ,a 4 ^2-- For = 1 we have the only term a 4 pM . The coefficients Vi could be straightforwardly 
obtained in terms of I, m, n, V , m' , n'. 

As a conclusion, we have fitted our results over the whole range of p 2 according to the following 
formula, where the Cj's are constants independent of p 2 : 



_2_[4] f n 2 ii A \\ 2 n 4 T)[6] 

Q(p 2 ,aV 4 UV 6, 1 -) = Q(p 2 ! ) 0)+c 1 ^ + C2 ^- + C3 ^ + C4 aV 41 , (31) 

\ P I P 

with indeed small x 2 S- We have also checked the validity of this expansion for the free propagator. 

The functional form used for Q(p 2 ,0, 0) does not influence significantly the resulting artifact 
coefficients. We can even avoid using any assumption about this functional form by taking all the 
values for Q(p 2 ,0, 0) as parameters which can be fitted 7 . 

This improved correction of hypercubic artefacts turned out to be particularly necessary for Z^. 
The results, including overlap-computed quantities, have already been presented in [38]. The raw 
lattice data for Z$ and Zy IOM exhibit dramatically the "half-fishbone" structure which is a symptom 
of strong hypercubic artefacts, and we recall that these effects are especially strong in the overlap 
case. After applying the pi 2 "] extrapolation method, eq. (|31j) . the curves are now perfectly smooth ; 
they do not either exhibit the oscillations which remain in previous methods. We will return later to 
the fact that Zy IOM is not at all a constant. 

Altogether we would like to recall the following hierarchy: first, the hypercubic artefacts are one 
order of magnitude larger for overlap quarks than for clover ones. Second, for both types of quarks 
the hypercubic artefacts for Z^ are one order of magnitude larger than those for g\. 

Let us stress that the discretisation artefacts we are discussing are all due to the QCD interaction. 
Indeed, our definition of Z 1 p{p 2 ) is such that it is equal to 1, as in the continuum limit, when interaction 
is switched off, and we have also g\ = 1 in that case. This illustrates that, in general, it may not be 
sufficient, by far, to extract the free case artefacts, or to use such prescriptions as replacing p^ by p M . 



4.2.3 Quantitative comparison with the "democratic" method 

We would like also to recall the quantitative comparison of our " pl 2n l extrapolation method", eq. 
piU with the more common "democratic selection" methods ; the latter method is carefully defined 
in [39]. This comparison is important, since almost all works up to now are using some variant of the 
democratic method, and since the difference with this method is crucial, as we show, to extract power 
corrections. 

Let us consider the overlap case. If we try to select, [39], the orbits which are in a cylinder 
around the diagonal with a radius 2tt/L, this is too restrictive anyway for our 16 4 overlap case, where 
we have only 22 orbits at the start. In order to have a less restrictive democratic criterion and to 
make a bridge with our own method, we will use the p[ 2n l's defined in eq. (|29j) . In our language, 
democracy can be translated as having a small enough ratio p^/{p 2 ) 2 . Momenta proportional to 
(1,1,1,1) and (1,0,0,0) have ratios 1/4 (minimum ratio, maximally democratic) and 1 (maximum, 
totally undemocratic) respectively. We then retain the "democratic" orbits defined by an intermediate 
p^/(p 2 ) 2 < .5. This leaves already only 7 orbits out of 22 for every j3. 

7 We have enough data for that. 
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In fig. [2 we plot for Zw, the result of this selection, as compared with our own method. Fig. [2 
clearly shows visible oscillations in the democratic curve demonstrating that the hypercubic artefacts 
have not been totally eliminated, while our treatment yields something perfectly regular. We prefer 
our own method for this reason and also because of the loss of information due to the rejection of 
"undemocratic" points, which leaves us with very few points. This appears crucial in the overlap case, 
where one is constrained to use small lattices. 

Then another very important aspect appears : while completely eliminating the hypercubic oscil- 
lations, we also considerably modify the mean value of the curve, as defined by an analytic fit ; for 
the present case, as seen on the figure, our curve is considerably lower, with a quite steeper descent 
than the democratic one. 



3.5 



2.5 



1.5 



(3 = 6.0 without hypercubic artefacts 

(5 = 6.4 without hypercubic artefacts 

p = 6.0 democratic selection 

p = 6.4 democratic selection 



I * J J, 



hi 



Figure 1: Comparison between the "democratic" method and our refined treatment for the 
overlap Z^(p 2 ) at 6.0 and 6.4. Above 2 GeV, precisely where one must study the power corrections, 
the difference is striking. 



4.2.4 The importance of optimising the elimination of hypercubic artefacts 

Let us stress then that, as is particularly visible from this comparison, the difference of the a 2 p[ 2n l 
method with the standard democratic method is not at all academic in the context of the study of 
power corrections and renormalisation constants. The difference with the previous versions of our 
method is also not negligible, as we have found. 

1) It is obvious that with the standard democratic method, we would obtain quite different 
results for power corrections and 0(4)-symmetric artefacts, and therefore for the resulting perturbative 
contribution. In fact, it has not even been really considered that could be affected by such power 
corrections. 2) Moreover, we observe that the power corrections, as well as the residual 0(4) symmetric 
discretization artefacts, extracted by the previous variants of our method for treating the hypercubic 
artefacts are not the same ; indeed, when using a previous cruder treatment for overlap action, we have 
found an important a 2 p 2 artefact, which disappears with the more refined treatment, and we were 
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also finding different power corrections (with a weaker condensate). The Wilson case shows similar 
spuriosities. This means that for a too crude treatment, some hypercubic artefacts can be spuriously 
mimicked as part of 0(4) symmetric discretization artefacts or continuum power corrections. This 
does not imply that the determination of power corrections is uncertain in this respect, but rather 
that it is very important to push hypercubic artefact elimination to the best to obtain the genuine 
continuum power corrections. 

Let us finally mention an interesting consequence ; as is well known, the values of Green functions 
at different momentum points in a Monte-Carlo lattice calculation are highly correlated, which should 
lead to very small x 2 /d.o.f. for fits describing the p dependence by smooth analytic expressions (by 
small, we mean well below one). It is not found so with too crude treatments of the hypercubic 
artefacts, because of the erratic oscillations which always remain in the latter methods mimicking 
statistical deviations. We observe however an impressive decrease of the x 2 down to its expected 
small value when we improve the treatment of the data, showing that we are now obtaining indeed 
very smooth functions as physically expected. 



5 Proof of the presence of "non canonical" artefacts : 



Problems for the determination of the continuum 



We now start, in the rest of the discussion, from the data obtained through the above treatment 
of the hypercubic artefacts. They still differ from the continuum by renormalisation and by 0(4) 
symmetric discretization artefacts. It happens that the determination of these artefacts is still harder 
in general than for hypercubic ones, (as for finite volume artefacts, which we estimate to be weak, 
see the corresponding short section below). 

Indeed, in practice, there is no similar unambiguous method to determine the 0(4) symmetric 
discretization artefacts. True, they manifest themselves by a residual variation with f3, and, in principle, 
we could study the variation with j3 at each momentum, and then extrapolate to the continuum. 
However, this requires too many momenta and too much accuracy, if we want really to extract the 
power correction from the extrapolation to the continuum. 

a) a general method for treating 0(4) scalar artefacts 

Then, a more practical method consists in assuming a prescribed analytical form for both the 
continuum and the artefacts, with some unknown parameters to be determined by x 2 adjustment. 

One has therefore to appeal to our a priori knowledge a) of the continuum, as function of p ; 
b) of the structure of 0(4) symmetric artefacts, as function of p and a, so that we could make fits 
with prescribed functions depending on a limited number of free parameters. For the continuum, 
this is exactly what is provided by the OPE, with the renormalisation constants and {A 2 ) as free 
parameters. For 0(4) symmetric artefacts, a standard idea is to recourse to lattice perturbation 
theory, in which the structure of artefacts is easily explicited. After all, this is what we have invoked 
for the hypercubic artefacts. As for 0(4) symmetric artefacts, the result is quite simple : in the case 
of a scalar function, and in the chiral limit, where there is no other dimensioned parameters than 
a and p, there can be no other artefacts than a n p n , with n > 0. Then we could work with a few 
parameters only. 

b) why it fails 




symmetric discretization artefacts. 



limit. 
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But now, there appears a very unfortunate circumstance. It is seen that this usual assumption of 
a perturbative structure of artefacts does not work at all, at least in general. We observe undoubtedly 
0(4) symmetric artefacts decreasing with p, i.e. for instance of the form l/p n and n > times some 
positive power of a. To show it, the best way is to consider the p dependence of Zy OM = gx/Z^, 
which should be momentum independent close to the continuum limit. Any p dependence is therefore 
to be attributed to artefacts. Now, the lattice data show quite clearly a very strong p dependence of 
Zy except at large p. This is true for clover action, and for overlap action as well, see Fig. Elfor 
overlap fermions. 




Figure 2: Overlap Zy OM as function of momentum, respectively for j3 = 6.0 and 6.4, and (3 = 6.6 
and 6.8, for the overlap action. The strong momentum dependence which is displayed must be 
entirely attributed to artefacts. 



The artefacts are similar ; they decrease monotoneously with increasing p for (3 = 6.4,6.6,6.8; 
they could be fitted by negative powers of the momentum squared ; for the overlap case, this is also 
true for /3 = 6.0 . 

c) consequences for determination of power corrections 

The necessity of taking into account such negative powers of p is very embarrassing, because, 
as we explained, we cannot distinguish the 0(4) symmetric discretization artefacts by the sole (3 
dependence ; we have to rely on the p dependence. Now, it is clear that we will have much difficulty in 
distinguishing the power corrections and the artefacts decreasing with p, and therefore to determinate 
the condensate. One of the problems is that when increasing the number of negative power terms in 
the description of the artefact, the tendency is to get alternating signs, and therefore rather unstable 
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results. There is no rationale as to where we should stop. In addition, we must consider that we 
have also as parameters in the fit the renormalisation constants Zb, which are practically chosen as an 
independent parameter at each (3. To add to the uncertainty, we observe that within the precision of 
the data, we can obtain equivalently good solutions by modifying the 0(4) artefacts with a correlative 
change in the Zb's. Then, in spite of many efforts, we have in fact not been able to extract stable and 
accurate values of (A 2 ) by this method, although we have got a clear signal that it is positive and 
sizable. Then, we have been able to fix the continuum power corrections only by exploiting particular 
circumstances, to which we devote the rest of the analysis, after a few words on finite volume artefacts. 

6 Finite volume artefacts 

Let us recall that in the Wilson case [38], we have found only very small volume artefacts, after a 
careful study with 8 4 ,16 4 ,24 4 lattices. Only the first points with smallest n 2 were showing some 
effect. 

We have not performed the same tests for the overlap action, because of the slowness of numerical 
calculation ; we have only ran on a 16 4 lattice. We are conscious that this is a possible weak point, 
since this volume is small, and we rely mainly on the analysis of overlap results for the determination 
of the condensate. So we think of extending this analysis to larger volume as soon as possible. 

We stress however that the consistency which is obtained for the continuum power corrections 
between the overlap results and the S-W (clover) ones with a greater volume 24 4 (see section EJ, is a 
further proof that the volume effects are not crucial for our purpose of determining the condensate. 

We must also observe that we can determine the condensate mainly from our large p data (p > 
4 — 5 GeV)(see 17. lj) . with large Lp, where volume artefacts are expected to be even smaller. 

6.1 Discretisation or volume artefacts ? 

The fact that the above 0(4) symmetric artefacts are observed at small p -and only at small p 
is deserving a special discussion, since it is so counter to usual expectation. Our ears are indeed 
accustomed to the dictum : "artefacts at small p, finite volume artefacts" , "artefacts at large p, finite 
spacing artefacts". Why are they not volume artefacts ? Since it is one important finding of our 
study, we collect here our arguments. 

Let us stress that, in the overlap case, the use of a small volume 16 4 at large (3, to which we are 
constrained, is not the reason for our finding of large 0(4) symmetric artefacts at small p. The first 
argument is that they are also seen in the clover case where we have tested in detail the smallness 
of volume artefacts; to reiterate, we have seen that only one or two of the smallest momenta seem 
affected by a volume dependence, while the artefacts we are discussing now are present over a large 
number of points . 

Moreover, let us stress that they have the typical behavior of discretisation artefacts, i.e. they 
decrease at large j3 and fixed p : indeed, Zy (which should be flat in the continuum limit) is flatter 
and flatter as beta increases at fixed number of sites 16 4 , i.e. as the physical volume decreases. More 
precisely, we observe that, as the volume decreases, we have more and more points where it is flat, 
i.e. where it has small artefacts : the respective number of points is at j3 = 6.0, 3 at j3 = 6.4,9 at 
j3 = 6.6, 14 at (5 = 6.8. This is exactly counter to a volume effect, for which the flatness should be 
obtained for n 2 larger that some fixed number. In other terms, for a volume effect, we would expect 
the artefacts to vanish at smaller and smaller values of p as j3 decreases (i.e. when the physical volume 
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increases). Instead, we find that they vanish around 4 — 5 GeV irrespective of j3 = 6.4,6.6,6.8 ; even 
more, at (3 = 6.0, with the largest volume, Zy never becomes flat . 

On the other hand, we have no understanding of why Zy IOM seems devoid of sizeable artefacts 
at large p : it must be regarded as an accident. But it must be underlined that this is not an accident 
specific to our particular problem. It is a well-known fact, on which any MOM practitioner is relying 
without being able to explain it : one measures the Z's in regions of large p, assuming that artefacts 
are small, although they would be expected to be large precisely there from lattice perturbation theory. 

7 Cases with small 0(4) symmetric artefacts with the 
overlap action at s = 0. Determination of the con- 
densate. 

As we have underlined as conclusion of section El paragraph c), the 0(4) symmetric artefacts we have 
found, when present, impede a clear determination of the power correction. It is fortunate that they 
seem absent or small in some cases. All correspond to overlap action ; then the results are convincingly 
consistent. Therefore, we concentrate on the overlap fermions in the present and following sections 
(|7I8|) . As concerns the clover action, we will not be able to have so compelling conclusions, but we 
show that they are at least compatible (see sectional). 



7.1 Large p (p ~ or > 5 GeV) 

At this point, we notice that the above difficulties may disappear first at large p, above roughly 5 GeV 
(therefore at (3 = 6.4,6.6,6.8), because Zy is constant in this region, especially for the overlap action 
; therefore it is suggestive that the artefacts of Z^ and g± are small there, barring for the unprobable 
eventuality that and g\ would happen to have exactly the same artefacts. We then fit Z^ and g\ 
with a formula containing no 0(4) artefact,, i.e. we take only the continuum expression times the 
renormalisation factor z^, at fixed /?, and we obtain then a very encouraging conclusion. Choosing the 
window of p within which Zy is well constant, we obtain at each j3 = 6.4,6.6,6.8, and for Z^ and 
gi, i.e. for six independent data, almost the same condensate value ; we quote a common fit (see the 
FigEJ) to the three Z^, with p m in = 5,5,4 GeV respectively, corresponding to the respective p where 
Zy is beginning to be flat : 



(A 2 ) ~ (3.1 ± 0.3) GeV 2 (32) 

large and positive, with a rather small error, and quite consistent with what is found in the gluon 
sector (see below). The x 2 /d-o.f. = 0.07 is very small as expected from very correlated data. Let 
us recall that this value of (A 2 ) is obtained with the convention that the Wilson coefficient of the 
operator is expressed in terms of a^f OM . 



A remarkable feature of this region of momentum is that both Z^ and g\, and therefore Zy too, 
are almost independent of f3. We take it as a further indication that artefacts are accidentally small 
there. 

Let us reinsist that the value of the condensate given in eq. (|32j) corresponds to the convention, 
always followed in this paper, that the Wilson coefficient of t he operator, calculated only at leading 
order, is expressed in terms of a MOM . The choice of a MS , would lead to an appreciably higher 
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Figure 3: Renormalised overlap 2L, at (3 = 6.4,6.6,6.8 and large p, compared with the purely 
perturbative result . The fit of the solid line is made with the perturbative part and only the 
condensate term , in the windows where the respective Zy OM become flat (solid line). One 
finds (A 2 ) = 3.1 ± 0.3. Note that is still far from being flat. It has a monotoneous decrease, 
of around 6% over the range, mainly due to the condensate. The dotted line is the purely 
perturbative result, clearly inconsistent with the data. Normalisation is made at 10 GeV. g\ is 
described by a similar fit. 



value (larger by around 70%). However, what is important is that the power correction by itself is 
well determined by our analysis, almost independently of such a change. Indeed, let us replace the 
OPE expression for the power correction by a simple power, without logarithms corresponding to a s , 
while maintaining the full perturbative expressions for the perturbative part. The fit then gives for the 
coefficient c of the power term c/p 2 : 

c = (0.767 ± 0.083) GeV 2 (33) 

when using a MOM in the perturbative part, and : 

c = (0.844 ± 0.083) GeV 2 (34) 

a small change indeed, of only 10%, reflecting the small change of the perturbative part, which is 
calculated by theory to a great accuracy. 



7.2 Z^p over the whole allowed range of p 

In addition, we observe another fact : Z$ - but not g\ - is strikingly independent of (5 over the full 
range of p. We have no explanation for that, but we can at least interpret this as meaning that Z$ 
is free of 0(4)-symmetric artefacts over all this range. And indeed, we can fit Z^ on the full range 
p > 2.6 GeV allowed for the perturbative calculation, and the four j3's with : 

(A 2 ) = (2.73 ± 0.21) GeV 2 (35) 
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See the figure |4] 
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Figure 4: Overlap Z^ at /3 = 6.0, 6.4, 6.6, 6.8, renormalised at 10 GeV. The fit for Z^ extends to 
the full allowed range of momentum, with (3 = 6.0 included, and a close value of the condensate 
{A 2 ) = 2.75 ± 0.2. The dotted line is the purely perturbative result. 



We thus obtain a remarkable similarity of the condensate with the previous value. We also obtain 
similar results by varying the window of the momenta, provided the lower limit is not pushed beyond 
5 GeV, and also by selecting various triplets of /? values. This seems to support the consistency of 
our assumptions about artefacts. 

If we want to still improve the agreement with the large p analysis, we can introduce a further 
l/p 4 term, accounting for the possibility that we may not be sufficiently asymptotic to have a good 
description with the l/p 2 condensate alone. For the full window, and the four (3's, we obtain a very 
good fit with : 



(A 2 ) = (3.2 ± 0.3) GeV 2 



and with a small subleading term: 



(-1.1 ±0.4) GeV A /p A 



(36) 



(37) 



with a minus sign which explains that the value of (A 2 ) is slightly higher than in the preceding fits. 
It is then closer to the large p fit. 



7 3 %MOM 

From our analysis, we can deduce values of Zy OM standardly defined, i.e. free from artefacts 



Z 



MOM, large p 
V 



nrMOM, large p 



(6.4) 
(6.6) 



1.798 ±0.012, 

1.776 ± 0.003, Z$ OM ' large p (6. 



1.756 ±0.011 



(38) 
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"large p" superscript being to recall that we have selected the value at large p, where we hope to 
have small 0(4) artefacts, in view of the observed flatness ; the precise value is obtained by a fit. For 
(3 = 6.0, we quote Zy IOM ' large p (6.0) = 1.878 ± 0.014, but Z^ OM is not yet flat at the highest 
momentum. Let us repeat that it is remarkable, and perhaps surprising, to observe such a constancy 
with j3, while one loop perturbation theory would predict a strong variation with f3. 

7.4 Comparison with lattice perturbation theory 

The fact that and Zy are very different from 1 (in fact not far from 2.) in the overlap case with 
s = may seem surprising. But in fact, already in lattice perturbation theory, the tendency is that the 
one-loop corrections to Z^ and Zy are large, because of a very large tadpole contribution to the self 
energy [29,31] (while g\ remains close to 1 as we find non pertubatively). The net effect is already large 
at (3 = 6.0, for the usual s = 0.4 : Zy = 1.247, and still larger for our s = : Zy = 1.444 ; we use 
Table 1 of ref. [29] for the analytical expressions (the definition of Z^ is different, but by a negligible 
amount); the numbers are quoted assuming a BPT boosted coupling with g 2 BPT = 1.68 (see ref. [31] 
under eq. (44)). Of course, what is surprising is first that the non perturbative determination (|38j) is 
still much larger, and, second, that it is almost independent of (3 over a large range. 

Note that the value found by ref. [30] for Za at (3 = 6.0, Za = 1.55, which should equate Zy, is 
also much larger than BPT ; however, this is not Z^ OM as measured directly ; it is the Za deduced 
from hadronic W-T identities. One may think of large artefacts which render different the results 
from various definitions of ZA,Zy ; our study below, section W2\ shows indeed the presence of such 
effects, very large at (3 = 6.0; but they are decreasing rapidly at larger (3, and our finding is that, at 
(3 = 6.8, Zy IOM is still much larger than the boosted lattice perturbation theory. 

8 Zy and Za- Consistency checks of the overlap action 
results. Complementary studies on chiral symme- 
try and artefacts 

To ascertain the soundness of our analysis, which may be surprising in several respects, we also perform 
several consistency checks of general nature in the overlap case ; indeed, some strong statements 
deriving from chiral symmetry can be formulated. 

8.1 Za/Zv in the non perturbative MOM scheme 

One should expect Z^ OM /Zy OM for the improved Green functions to be exactly 1, i.e. without 
any artefact in the chiral limit, and in the perturbative regime, i.e at large momentum. This derives 
from the exact chiral symmetry of the action, discovered by Luescher, combined with the choice of the 
improved current, or equivalently, of improved Green functions. However, one has also to assume the 
symmetry of the vacuum state, therefore the absence of spontaneous symmetry breaking. Therefore 
this result only holds at large momenta. We find indeed this to a high accuracy, fig. El Above 2 GeV, 
for the four /3's, we have 

Z MOM/ Z MOM = 1 tQ 

a very high accuracy. On the other hand, we see 
that for lower momentum i) the ratio differs from 1, and ii) it depends on j3. The change of regime 
is rather abrupt. Observation i) has also been made for domain wall fermions at j3 = 6.0 [41]. 
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for overlap fermions and small quark mass at the four (3 = 
The ratio is 1 above 2 GeV, but differs from 1 and depends on (3 below 2 GeV . 



8.2 Za from a hadronic W-T identity against Zy I0M 

In principle, once artefacts and power corrections have been eliminated, the various ways of defining 
the renormalisation constants can be related by perturbation theory. In the case of Zy or similar 
cases, for the same action, implying identical finite parts, they should even be equal. Moreover, Za 
should be equal to Zy from the chiral symmetry of the overlap action. 

We verify this statement, as a highly non trivial consistency check of our treatment, by measuring 
Za from a standard hadronic Ward-Takayashi identity. Z^ 1 = m q /p, where : 

p = 1/2 < d A , P 5 > I < P 5 , P 5 > (39) 

In contrast to Zy OM,lar9e p , Z^ 1 presents a very strong variation with (3 : at 6.0, the situation seems 
hopeless, with Z^ 1 ~ 3, but the decrease is rapid : 

Zjf 7 (6.0) = 3.03, Z^\<6A) = 2.046, Z^(6.6) = 1.90, Zf 1 (6.8) = 1.78 (40) 

and one reaches finally a level close to Zy IOM ' lar9e p . Moreover, Z^ 1 (a) is linear in a 2 to a very 
good precision : 

ZY\a) = 1.65 + 5.52a 2 , (41) 

which shows that the difference Z^ 1 — Zy OM,lar9e p is indeed a discretisation artefact, as it should, 
and one of the most canonical species, since we expect precisely chiral symmetry breaking to be at 
most 0(a 2 ). Indeed, on hadron states, the Ward identity is valid up to 0(a 2 ), and Green functions 
at large p, where we measure Zy IOM,lar9e p , have also only 0(a 2 ) artefacts from chiral symmetry 
arguments, in the chiral limit. 
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9 Check through the Wilson clover 



In view of the rather paradoxical situation which allows us to determine (A 2 ) from the s = overlap 
action, it is important to reascertain this analysis by a study of the standard clover action, to check 
whether we can obtain a similar continuum result. Let us stress indeed that from renormalisability, 
the whole continuum function must be the same up to an overall constant for all the versions of the 
action. A priori, as we have said, the Wilson case may seem hopeless because there are large 0(4) 
artefacts even in ; and even at large p, Zy is not so flat 

However, one can benefit from the knowledge gained in the overlap case. In the end, the situation 
appears in fact somewhat similar for the clover action as concerns 0(4) symmetric artefacts: the main 
burden of non canonical artefacts appears to concern g\, and we can obtain a good description of Z^ 
with a minimal a?p 2 canonical artefact, and just one a/p 2 term. Indeed, we obtain : 

(A 2 ) = (2.4 ± 0.3) GeV 2 (42) 

with artefacts 

- (0.005 ± 0.0015) a 2 p 2 - (1.9 ± 0.4) GeV 3 a/p 2 (43) 

In the present case, the role of these artefacts is crucial to obtain the condensate. In fact, before 
extraction of these artefacts, the clover data show a rather flat behavior at large p, or even an increase 
at 6.8. The term with a/p 2 considerably improves the fit, with a \ 2 divided by 3, and brings in a 
value of the condensate quite close to the one obtained with overlap fermions . These results seem a 
signal that we have obtained a rather accurate treatment of artefacts. Of course, one may be worried 
of having introduced a non canonical artefact, which was not necessary in the overlap case ; yet, we 
must remember that there is no logical reason why such terms should be absent (in the overlap case, 
they are present anyway in gi). 
With a subleading term 

oc 1/p 4 (44) 
we obtain a still somewhat better agreement with the overlap condensate : 

(A 2 ) = (2.83 ± 0.35) GeV 2 (45) 
The subleading term is also of same sign and comparable magnitude as for overlap action : 

(-1.85 ± 0.85) GeV A /p A (46) 

Note that the consistency of the Wilson 24 4 results with the overlap analysis is comforting the 
idea that the volume effects are not too important in the overlap case, although certainly they are less 
compelling because of the need to introduce the a/p 2 artefact. 

10 Consistency with the gluon analysis 

To stress the overall consistency of the analysis, one has also to consider the agreement with the 
previous gluon analysis. We quote the result of a combined fit of (A 2 ) and Aqcd to the gluon 
propagator and the symmetric three-gluon vertex [19] : 

(A 2 ) = (3.6 ± 1.2) GeV 2 ( a M OM ), (2.4 ± 0.6) GeV 2 ( gluon propagator ).... (47) 
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with Aqcd = 0.233 ± 0.028 GeV, and three-loop anomalous dimensions (the use of our present 
Aqcd = 0.237 would lower somewhat the condensate values). The fit has been done with the same 
convention as for quarks, that the Wilson coefficient (calculated only to one loop) is expressed in terms 
of a^ OM . (We could get free from this convention by comparing the magnitude of power corrections 
themselves as in eq. The coincidence of the values of (A 2 ) with the quark value is highly 

significant since it concerns the continuum function, extracted by a series of various, independent, 
manipulations committed on the gluon and quark Green functions, and since these continuum functions 
are related only thanks to the OPE. 

The quark measurement has much smaller statistical errors. This is simply due to the fact that 
in the gluon case, we have left free the value of Aqcd, and moreover, the value of the condensate 
depends strongly on this value, whence the large errors in the gluon case. On the contrary, in the 
quark case, we have chosen a fixed Aqcd = 0.237 GeV. Indeed, in this case, the dependence on 
Aqcd is rather weak. Then it is useless to try to determine Aqcd from the fit, and on the other hand 
the value obtained for (A 2 ) remains well determined if we allow some variation in Aqcd- Another 
advantage of quarks is that we can reach the accuracy of four loops in the theoretical expression for 
the purely perturbative part. In the gluon case, such an accuracy is only possible for the asymmetric 
vertex, but in that case, the leading term in OPE is not given by (A 2 ) [20]. 

11 Conclusions and discussions 

11.1 Physical results and evidences of various artefacts 

Let us summarize our results in the following few points : 

• First of all, and this is the main point, let us stress that we have finally obtained a rather 
non trivial confirmation of the validity of the OPE in the non gauge-invariant sector of lattice 
QCD (treated numerically). The virtue of OPE is that one can describe the departure of all 
the various Green functions from the perturbative approximation at large momenta with the 
same set of expectation values. We have now obtained consistency not only between gluon 
Green functions, but also with the quark sector. This is highly non trivial. Let us recall that 
there is often some doubt raised about OPE itself, and the possibility is sometimes considered 
that power corrections could be present not corresponding to an operator v.e.v. The final 
consistency of the determination of the condensate from three distinct type of Green functions 
strongly suggests that this is not the case, at least for leading power 1/p 2 corrections. We 
do identify an OPE power correction, consistently related to A 2 . (A 2 ) is found consistently 
large and positive. The precise magnitude of (A 2 ) is affected by an important uncertainty, due 
to the low accuracy of the theoretical calculation of the Wilson coefficient. But the power 
correction (i.e. the product of the coefficient and the condensate) is well determined by the 
lattice analysis, and the ratio of the power corrections in the various Green functions is actually 
as expected from lowest order OPE. 

• It turns out that the lattice discretisation artefacts are unusually sizable in the quark propagator 
Zip, but a clearcut distinction must be made between hypercubic artefacts, which are gigantic, 
but can be efficiently eliminated, and the 0(4) symmetric ones, which are not so catastrophic, 
but that we have not been able to handle systematically. 

• We believe that we have been really efficient in getting rid of the hypercubic artefacts thanks 
to our (improved) method of "restoration of 0(4) symmetry". 



22 



Once these artefacts have been subtracted, the overlap Zy(p 2 ) = Z^{p 2 ) / 'g\{p 2 ), which should 
be independent of p except for artefacts, is very close to a constant at large p > 5 GeV. This is 
far from trivial and supports the statement that we have no remaining 0(4) symmetric artefacts 
in this specific region. This is directly supported by the near constancy of the quantities as 
function of (3. 

Moreover, also in the overlap case, but for Z 1 p{p 2 ) only, the same statement of constancy with 
(3 extends down to the lowest momenta ; this leads to suspect that 0(4) symmetric artefacts 
are small in this case over the whole range of p. We are unable to explain these two special 
situations . 

Let us recall that a certain flatness Zy(p 2 ) at large p is also observed with the Wilson quark 
action, although it is not so good. Let us recall also that this region is the basis for the standard 
determinations of MOM renormalisation constants, usually with the Wilson quark action (see 
for example [42], which presumes that discretisation artefacts are not large there. Let us then 
recall that we have no theoretical argument supporting the statement that they are not large. 
Quite the contrary. If any, the theoretical arguments would suggest them to increase with p. 
The support is purely empirical. This is embarrassing if we aim at precision determinations. 

Zy of the overlap action at s = is large, around 1.8, in qualitative agreement with BPT per- 
turbation theory which finds large self-energy contributions in Z^(p). But it is still much larger 
than the expectation, and the lack of dependence on (3 is not understood from perturbation 
theory. 

Considering the cases where the 0(4) — symmetric artefact-free results are supected to be 
small, we try to fit them by OPE, i.e. by the four-loops perturbative contribution plus the 
< A 2 > condensate contribution computed to leading logarithm. The overlap Z^ and g\ at 
large p > 5 GeV allow for a good fit for (3 = 6.4,6.6,6.8, leading to a consistent < A 2 > not 
far from 3 GeV 2 . The overlap Z,^ also allow for a good fit for the whole range p > 2.6 GeV , 
including in addition (3 = 6.0. The < A 2 > condensate is consistent with the former value. A 
very small l/p 4 term still improves the consistency. 

In the other cases, namely in g\ and Zy for p lower than p ~ 5 GeV , the 0(4) — symmetric 
artefacts become large, especially at small p, and in fact, they increase regularly from large to 
small p. This trend, which is also contrary to the expectation of lattice perturbation theory, 
clearly indicates a non perturbative origin. In fact, one important conclusion of our study is 
the existence of these very large non perturbative artefacts at small p due to discretisation. 
These are very embarrassing for any analysis of the Green functions, as we comment below. 

Let us repeat our strong conviction that these low p effects(under 2 GeV) are indeed discreti- 
sation artefacts and not volume artefacts. They behave quite counter to volume effects, as we 
have extensively argued. 

A short study of Z^ OM /Zy OM shows that similar artefacts are still present in a ratio where 
they would be expected to cancel if one applies naively the exact chiral symmetry of the lattice 
action. They are clearly seen as being discretisation artefacts because they are reduced at larger 
(3. They seem to be present on top of an actual continuum effect, small but visible, which 
could be due to continuum chiral symmetry spontaneous breaking. One can suspect that the 
artefacts themselves are connected with the spontaneous breaking of chiral symmetry. 

It seems logical that such chiral-symmetry violating artefacts, as well as the continuum effect, 
be only forbidden at large p, if they are connected with spontaneous breaking effects. Indeed, 
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it is there and only there that such effects fade away. Therefore, it is consistent with this 
interpretation that we find Z^ OM /Zy OM = 1 at large p to a high precision, without (3 
dependence. 

• Of course, the clover (SW) action has the great advantage that it does not present the same 
very constraining limit in volume as the overlap one ; however, it leads to less compelling results 
than the overlap one, because we have not found here the same particular situations where 0(4) 
artefacts can be neglected ; the small p artefacts, which render so difficult the determination 
of power corrections, are present in and not only in g\. As explained above, by including 
more and more terms to describe these artefacts, we destabilize the numerical value of (A 2 ). 
Stopping with the first term a/p 2 , we obtain consistency with the overlap results. 

• As to comparison with other works on the quark propagator, the question of the presence and 
magnitude of power corrections is a crucial test of the precision obtained in the treatment of 
Green functions : the condensate value (A 2 ) should be independent of the action and of f3 with 
due renormalisation. In our opinion, safe and accurate extraction of power corrections requires 
a very large range of momenta, and therefore a large range of (3 8 ; indeed, to use a rather large 
range of momenta at a fixed j3 would be dangerous because of the periodicity of the lattice ; 
the large p behavior would then be highly dependent on empirical redefinitions of the momenta, 
aiming to remove empirically lattice artefacts at very large p ~ 10 GeV . Of course, a crucial 
question is whether one can work with a j3 as large as 6.8 with our lattice size 16 4 - the possible 
size is indeed strongly restricted for the overlap action. It is our conviction, for reasons which 
have been explained in detail. Another concern is that, according to our experience, to extract 
the real power corrections, one needs a particularly careful elimination of hypercubic artefacts. 

• The resulting value of (A 2 ) from the OPE analysis of lattice QCD data should be compared 
to tentative estimates made by various authors within analytical approaches. One will find 
abundant references in the paper of Dudal et al., ref. [43]. It is clear that this comparison must 
take care of the precise definition of the condensate, as regards for instance renormalisation. 

11.2 Systematic errors on (A 2 ) 

Of course, we are making many assumptions which introduce uncertainty in the value of (A 2 ). Recall 
that we do not claim to determinate it only from the present study, since we have several previous 
determinations from the gluonic sector. So we can also appreciate systematic errors from the consis- 
tency we obtain with these previous estimates (see above, section ITU]). In fact, as we have observed 
in section [^1 the gluon propagator determination and the one from a s are notably different, but 
the values are affected by very large errors, and are compatible with the present ones. From all the 
determinations, we could conclude that the systematic errors do no seem to exceed 1 GeV 2 . But 
there is in fact an important source of systematic error, which is explained below, and which cannot 
be estimated by comparison with gluons because it is present in both: it is the fact that the Wilson 
coefficient is calculated only at low order. Then, it remains useful to discuss the sources of errors 
inside the quark sector itself, for which anyway the conditions are intrinsically very favourable. 

8 The lattice group in Adelaide has performed extensive studies of the quark propagator (see references 
above). Recently, the same group has extended its analysis to overlap action [11] and a special gauge 
action, with several /3's, corresponding to a -1 ~ 1 GeV^ 1 — 2 GeV -1 , yet notably lower than our largest 
cutoffs. 
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Since we do not claim to do phenomenology, but rather an exercise in quenched QCD, we have 
not to bother about the quenched approximation, which is also supposed for gluons. Chiral limit is 
assumed on the theoretical side, for instance to calculate the Wilson coefficents. Now, of course, we do 
not work at zero quark mass on the lattice. We have not tried to do a systematic chiral extrapolation 
on the lattice data, which would only lead to increase the statistical errors. We observe that Z$ seems 
very weakly dependent on our set of masses, which means that this limit is not a priori a problem 
at the smallest mass (at which we have made all our OPE analysis). Anyway, we can discard any 
catastrophic effect at very low quark masses through the consistency with the quenched gluon data. 

Let us now pass to more relevant effects. 

Some come from the treatment of artefacts, for which we lack of theoretical basis, some are 
relative to our description of the continuum, which, although based on a much stronger theoretical 
basis, involves necessarily approximations. 

We do not return to finite volume artefacts, since we have nothing quantitative to say about their 
magnitude in the overlap case. For what concerns hypercubic artefacts, we may have an idea on the 
error remaining in their treatment by the variation observed with two variants, with a slightly different 
description of the continuum limit : 

(A 2 ) = (3.23 - 3.0) GeV 2 (48) 

Whenever large 0(4) symmetric artefacts are present, we are compelled, as we have seen, to rather 
arbitrary assumptions on their structure, since we have concluded that we cannot rely on lattice 
perturbation theory. The correlated errors in our determination of the "non perturbative" artefacts 
and of the condensate seem very large, as judged from the range of values obtained in various fits, 
to such a point that we have renounced to extract any number in this case. Therefore, we consider 
only the case where we have strong hints that the artefacts are small, in which case we have mainly 
to consider uncertainties in the continuum description. They are themselves of two origins : the 
perturbative calculation of Wilson coefficients ; the non perturbative aspect, i.e. the relevance of the 
OPE expansion, the enumeration of operators... As to uncertainties in perturbative calculations : 

i) We have checked that computing the perturbative contribution to third or fourth order in 
perturbation does not change significantly the estimated condensate (only 7% of change). Another 
test is to reexpress the series in terms of ctj^g instead of ckmom- Zip changes by less than 1% with 
various prescriptions. We can thus assume that the perturbative contribution has been expanded far 
enough. 

ii) As we have explained in subsection 12.31 the problem is much more important for the Wilson 
coefficient of the A 2 operator which has, on the contrary, only been computed to leading logarithm. 
A sign of this problem is seen by changing «mom into ojjg. A change of «mom into reduces it 
by fourty percent a s (10GeV), and more for the smaller momenta; whence a reduction of the Wilson 
coefficient by 50% in average. Through a conspiracy with the smaller change in the perturbative 
contribution, this change amounts to an increase of the resulting condensate by 70%. Of course, 
similar effects are present for gluons, and, as we have explained, the ratio of condensates obtained 
from quarks and gluons will remain the same. More importantly, one must be aware that the power 
correction by itself remains well determined; what is not well determined is the translation of the 
power correction into a (A 2 ) condensate value. Indeed, this translation depends on the theoretical 
evaluation of the Wilson coefficient, which is not accurate at present. 

As to the properly non perturbative aspect, we may think of two sources of uncertainties in 
determining (A 2 ). The one stems from other operators which could enter with the same power in 
the OPE expansion. However, we do not find any such operator contribution in in the chiral 
limit. We could have some contamination since we are not exactly in the chiral limit, but it must 
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be small, since we observe a very weak quark mass dependence. The other could be the possibility 
that p is not sufficiently large for the leading correction 1/p 2 to completely dominate over next ones. 
This possibility is represented by the 1/p 4 term and the fact that it is small but non zero in the fits 
shows indeed that we are not completely asymptotic at such large momenta ; although it is very small, 
around 1% at 3 GeV, it leads to a change of around 15% in (A 2 ). On the other hand, with this 
term included, we find a very good stability of (A 2 ) over the large range 2.6 GeV < p < 10 GeV, 
when varying the fitting window, which suggests that we have correctly accounted for the small 
subasymptotic effects. Actually, this term could also mimick a neglected logarithmic dependence ; 
indeed, passing from a MOM to a MS as expansion parameter as explained above, the 1/p 4 term passes 
from — 1.1/p 4 GeV 4 to — (0.65±0.34) GeV 4 /p 4 , therefore there is an appreciable variation, although 
the sign and order of magnitude are encouragingly stable. 

On the whole, the dependence of the Wilson coefficient on the scheme for a s seems the most 
worrying source of uncertainty, yet it can be solved soon. 

Another concern is the value used for Aqcd- Let us vary by ±10% our = 0.237 GeV. We 
find : 

= 0.215 - 0.260 GeV -► (A 2 ) = 3.45 - 3.02 GeV 2 , (49) 

a quite moderate change indeed. There is naturally an increase for decreasing Aqcd because the 
larger power correction compensates for the slower falloff of the perturbative part. 

11.3 General consequences for lattice studies induced by the ob- 
served discretisation artefacts and power corrections. Ac- 
curacy on Green functions. Renormalisation in the MOM 
approach 

True, the direct object of this study has been to verify the consistency of our OPE analysis of lattice 
data by extending it to the quark sector, and then comparing with the previous analysis of the 
gluon sector; and thereby to assess the soundness of our statement of large power corrections in 
"elementary" Green functions as well as of their interpretation in terms of the non gauge-invariant 
condensate < A 2 >. 

However, one must be aware of the strong consequences of this study, as well as of the preceding 
ones, on general problems, especially in precision studies : 

-1) Power corrections. First, we have a problem independent of the discretisation of the action 
: the presence of the power corrections has the effect of modifying the estimate of quantities defined 
in the perturbative regime of QCD, when one attempts to extract them from the numerical measures 
done on the lattice. Indeed, such power corrections, of non perturbative origin, must be necessarily 
subtracted from the Green functions to get the perturbative contribution, the only one which presents 
a universal character since it can be translated from one renormalisation scheme to another. 

In the gluonic sector, the perturbative contribution leads to a determination of Aqcd the 
necessary extraction of the 1/p 2 power correction induces a striking modification in the value of 
Aqcd. as was found some years ago [17]. Of course, this very large change of Aqcd corresponds to 
much more moderate corrections on the Green functions themselves, yet they may amount to several 
percents (10% at low points). 

Now, let us recall that for renormalisation constants of the quark sector, one is also most often 
looking for the ones defined in the perturbative regime. Indeed, only such perturbative renormalisation 
schemes can be connected between one another by analytical calculations, and also connected to 



26 



Wilson coefficients in order to produce physical quantities. But we predict from OPE, and we have 
indeed observed in this article, that power corrections of the same magnitude affect the quark sector, 
in particular and various vertex functions. From the Ward identity, it happens that for quantities 
like Zy, the A 2 power corrections cancel between the vertex function and Z^. But this is not true 
in general : in Z$ and Zp for instance, this cancellation does not occur. In such cases, the power 
corrections must be subtracted, in principle. Of course, one may wonder whether this is practically 
important. It depends on the accuracy we want to obtain. If we aim at a precision of a few percent, 
certainly we do require to take them into account, since they reach several percent around 4 GeV , 5% 
on Z^p in the Wilson case. Now, one often claims to go below 10% with dynamical quarks ; then, such 
effects are deserving of consideration. 

What simplifies somewhat the problem raised by power corrections is that not only they are 
independent of the chosen discretisation of the action, but they are often related to the same (A 2 ) 
condensate, at least as regards the dominant power. Once (A 2 ) has been confidently determined by 
one analysis, it can be used in others, the respective contribution to the various Green functions being 
obtained through the lowest order Wilson coefficients of the OPE. 

-2) the 0(4)-invariant discretisation artefacts that we have found set a more general difficulty, 
and one which is more embarrassing, in view of our lack of theoretical control. We have no reason 
to suppose that the special case we have studied is especially catastrophic. Yet, it shows already an 
embarrassing situation. 

a-The first step is to have control on hypercubic artefacts to a good accuracy. It is deserving 
of mention that this accuracy cannot be obtained by the standard method of selecting democratic 
points. Moreover, the simplest versions of our alternative method of " restoration of 0(4) symmetry" 
have not allowed to get a good accuracy. We have had to go further. This requires already a good 
deal of work. But, finally, it seems that we are able to produce a systematic and accurate procedure. 

/3-On the other hand, we have not obtained a safe a systematic general method to extract the 
0(4)-invariant artefacts while we have a clear proof that they are large at small p. This is disastrous for 
the extraction of the power corrections. This fact seems to be a problem specific to the quark sector 
: we had not found similar evidence for gluonic Green functions, although we cannot exclude totally 
their presence. We have not any theoretical mastering about their magnitude. It is our impression 
that the perturbation theory is of no help, since they happen in fact to be large at small p, and small 
at large p, counter to the expectation of perturbation theory. There is no explanation for the relative 
flatness of Zy we observe at large p, and which is the basis of most determinations of renormalisation 
constants. Note that this lack of explanation is a problem not only for us, but for all the MOM 
determinations of renormalisation constants. 

One may wonder whether the choice of the quark action may help. The use of overlap action 
seems to introduce larger hypercubic artefacts, but this is not a decisive obstacle, as we have explained. 
On the other hand, the fact that in certain cases it presents small 0(4)-invariant artefacts seems an 
important advantage. Yet, it is weakened by our lack of understanding of the underlying reasons why 
it is so. Calculations with other values of p = 1 + s may reveal instructive. 

Taking into account these uncertainties coming from 0(4)-invariant artefacts, which have appeared 
negligible only for Z^ and only in the overlap case, but not at all for Zy, it is possible that the MOM 
scheme may reveal not very practicable for precision calculations, although it is appealing by its 
simplicity, and quite efficient for ordinary purposes. Of course, at this point, one may think of the 
method of the ALPHA collaboration as a complementary one, technically difficult, but which allows 
a very clean treatment of discretisation artefacts by using on-shell quantities, and also allows to work 
at very high energy scale, therefore eliminating power corrections (see for example, for Zy, [44]). 
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A Perturbative expansion to four loops in the full 
MOM scheme 

Our aim is to express Z^(p 2 ), in terms of a^ OM , defined by the triple gluon vertex at symmetric 
momenta. More precisely, we are looking for the renormalisation group improved expression, which 
resums the large logs log(p 2 //U 2 ) at large p. The expression then takes the form of a series in 
a I f OM {p) with pure number coefficients. This series can be obtained from the knowledge of the 
anomalous dimension of in terms of a^ OM , and of the (3 function of the MOM scheme. This is 
possible in the Landau gauge to the order of four loops included, thanks to the papers of Chetyrkin 
and collaborators, ref. [25] and [26]. 

First, one has in section 4.2 of the first paper the expansion of 9ln ^^i ^ as a series in a s of the 

MS scheme to four loops (in fact the authors consider the inverse Z 2 (/U 2 ) = (Z^,(^ 2 )) -1 , so we have 
to invert their formula). Since there is no a s term in the Landau gauge, to reexpress the series in terms 
of the a s of the MOM scheme at four loops requires the expansion of <Jf S in terms of a^ OM only 
to order three included. This expansion is provided by the section 5 of the second paper, by inverting 
the first equation of this section. On the other hand, we need the MOM (3 function. It is given in 
the second paper, section 6, at three loops, which is also sufficient to calculate the renormalisation 
improved series for Z^(p 2 ) at four loops included ; indeed, the four loop coefficient of the (3 function 
enters only as a factor of the one loop anomalous dimension of the fermion 70 , which is zero in the 
Landau gauge. We then obtain the following expansions : 
-Z^ is expressed in general as : 

^ L , Poll ~ 2A70 



Zl ert = j | +n 



4nf3 2 



n - 2 f (A)7i-2A7o) 2 ^ 8/3 2 7o - Pofolo - Wi7i + 2/3 2 72 
+ 0.5a — — ■ „ h 
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+ { ~ 16 « + 24/3 /3 1 2 7 2 (-2/3 1 + 7l ) 

+ 2/3 2 /3i7o(-16/3 2 + 3/3 270 + 24ft 71 - 6 7l 2 ) 
- 2/^(2/3370 + foil + 4/3i72 - 37172) 



+ f3l (l6/3 2 7 i - 3/3 27o7l + 27? + 4/3i (2/3 27o - 3( 7 2 + 7072))) + 4/3 5 73 ) J (50) 
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-For rtf = in the Landau gauge, we get : 

7o = 0.0 7i = — 72 = -94.7943 

73 = 14503.7 (51) 

-a , which represents here q;mom(p) (defined through the three-gluon vertex at symmetric mo- 
menta), is given, under the same conditions, by : 

■e> - h - 8 4(^ + w? ( 2 ;l + ie 4 ((lnt)2 - tat - '>) (52) 

+ 7^1 f-32vr^(lnt) 3 + 80vr^(lnt) 2 



(/WV PI $> 
12* fofhlh - 647r/3f 27r/? 2 /? 3 - 16tt/3 3 

$ 11 /?o 3 

where t = ln(p 2 /A 2 AOM ) (in the MOM case, one stops at the 1/i 3 terms, and where 

(3 = 11 fj l= 51 /3 2 = 3072 (10%) (53) 

B A 2 Wilson coefficient in the quark propagator 

In order to renormalise the bare quark propagator in Eq. (|15j) we will define the following two 
renormalisation constants, both in the momentum subtraction (MOM) general scheme : 

V -ip' + mjp 2 ) / ^ ' 

<W = o ; , ox (54) 



^V)/ p2 + mip2) V Spe'nW / 



The constant in the top of the l.h.s of Eq. (|54j) includes the non-perturbative contributions to the quark 
propagator to let it take the tree-level value all over the energy range (not only in the perturbative 
regime). We renormalise Eq. (|15|) by multiplying by that of the bottom (this purely perturbative 
MOM renormalisation constant, computed to four loops in ref. [25], is presented in the appendix lA|) 
; as to the Wilson coefficient of A 2 , we calculate it at leading RG order. 




S pe M+rt - 2 4(N 2 -l) 6a > b + 



ZTV 1 ) , d {£>°M)<(A*) K M> \ 

where (A 2 )r([i 2 ) = Z~^{n 2 ) : A 2 ^ : 9 , and where we only write the leading terms in j>. Concerning 

d (jp> Qt (/ jt )) • tne same procedure used for gluon Green functions in refs. [18-20] is in order here. 

9 the : • • • fixes the normal order substracting the additive divergencies in A^ are (see discussion in 
ref. [40]) 
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Then, from Eq. (|55|), multiplying by Z^ eTt (fi 2 ) , and taking logarithm derivatives on fi in both sides, 
we obtain the following RG equation: 

^ wr M + ^ d ^^ (56) 

that is identically satisfied with : 

d ($ M ^) = d{ha{p)) (^y) "° ' (57) 

where d((l, a(p))) is in fact beginning with a(p) and 7^2 is defined through : 

and where all the involved one-loop coefficients are well known: 

7$ = ^' 70 = 0, /% = 11. (59) 

On the other hand, d (1, a(p)) can be obtained by computing the only diagram involved in the "OPE 
business" for this case, 



d(l,a(p))= {N ^ l) g 2 (p) (60) 



^ ) ^ 



Then, by applying the results from Eqs. (|57B59j) to Eq. (|55|). we will obtain the final result Eq. (fTo^ 
to be used for the fits. 
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